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Nonlinear instabilities of multi-site breathers in 
Klein-Gordon lattices 


By Jesus Cuevas-Maraver, Panayotis G. Kevrekidis and Dmitry E. 

Pelinovsky 


In the present work, we explore the possibility of excited breather states 
in a nonlinear Klein-Gordon lattice to become nonlinearly unstable, even 
if they are found to be spectrally stable. The mechanism for this fun¬ 
damentally nonlinear instability is through the resonance with the wave 
continuum of a multiple of an internal mode eigenfrequency in the lin¬ 
earization of excited breather states. For the nonlinear instability, the 
internal mode must have its Krein signature opposite to that of the wave 
continuum. This mechanism is not only theoretically proposed, but also 
numerically corroborated through two concrete examples of the Klein- 
Gordon lattice with a soft (Morse) and a hard (cp 1 * * 4 ) potential. Compared 
to the case of the nonlinear Schrodinger lattice, the Krein signature of 
the internal mode relative to that of the wave continuum may change 
depending on the period of the excited breather state. For the periods for 
which the Krein signatures of the internal mode and the wave cont inuum 
coincide, excited breather states are observed to be nonlinearly s table. 


1. Introduction 

The study of anharmonic modes constructed out of a few excited lat¬ 
tices sites in nonlinear lattice dynamical systems is a broad and diverse 
theme of research that emerged in the physics literature through the work 
of [211131]. Subsequently, the mathematical proof of the existence of such 
modes in [20] not only placed such states on a rigorous mathematical foot¬ 
ing, but also gave a deep set of insights towards their potentially generic 
nature. This, in turn, led to a considerable growth and diversification of 
interest into these modes over the last two decades, eventually culminat¬ 
ing in a wide array of reviews on both the methods of analysis of these 
modes, as well as on their diverse applications PEUnUIH]. Noting only 
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some among the many areas where these modes have been impacting, we 
cite coupled waveguide arrays and photorefractive crystals in the realm 
of nonlinear optics mm, the denaturation dynamics of the DNA double 
strand in biophysics EDI. breather formation in eletrical lattices and in 
micromechanical cantilever arrays [33], Bose-Einstein condensates in op¬ 
tical lattices in atomic physics [23], as well as bright and dark breathers 
in granular crystals mm- 

Orbital and asymptotic stability of the fundamental (single-site) 
breathers were established by Bambusi mm- Spectral stability of more 
complicated multi-breathers were classified in the recent works [2U3TU29], 
depending on the phase difference in the nonlinear oscillations between 
different sites of the lattice (and the nature of the potential). The sim¬ 
plest spectrally stable multi-breather configuration includes a two-site 
breather. If the two sites excited in the anti-continuum limit are adjacent 
in the lattice, then the spectrally stable breather is in-phase (anti-phase) 
for the hard (soft) potential V [ TT1I29] . 

The main question we would like to consider is if the spectrally stable 
two-site breather is also stable in the nonlinear dynamics of the discrete 
Klein-Gordon (KG) equation ((D)- This is a part of the more general 
and broadly important question about whether spectrally stable excited 
(non-fundamental) states of a Hamiltonian system are nonlinearly stable 
or not. In a similar context of the discrete nonlinear Schrodinger (NLS) 
equation, the two-site breathers near the anti-continuum limit are not 
orbitally stable because the linearized spectrum admits an internal mode 
of negative energy (negative Krein signature) [26]. This internal mode of 
negative energy destabilizes nonlinear dynamics of the two-site breathers 
if a multiple harmonic of the internal mode eigenfrequency occurs in the 
frequency spectrum of the wave continuum that bears positive energy 

hd ag¬ 
in the context of the discrete KG equation (pQ), a similar concept of 
Krein signature can be introduced for the internal modes in the lineariza¬ 
tion of multi-site breather solutions [3 ED [ 22 ]. The discrete NLS equation 
appears to be a valid approximation for the small-amplitude weakly cou¬ 
pled multi-site breathers [25] , Therefore, it is natural to expect that the 
mechanism of nonlinear instability based on the resonance between mul¬ 
tiple harmonics of the internal mode eigenfrequency and the frequency 
spectrum of the wave continuum may also be observed for the two-site 
breathers near the anti-continuum limit. This conjecture is fully con¬ 
firmed in this work. However, depending on the period of the limiting 
breather, one can find parameter configurations, where the nonlinear in¬ 
stability can be avoided, because the Krein signatures of both the internal 
mode and the wave continuum coincide. In such situations, our analytical 
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and numerical results indicate that the two-site breathers are stable in the 
nonlinear dynamics of the discrete KG equation ([Tj). 

The presentation of our results is structured as follows. In section 2, 
we present the mathematical formalism of the problem. Upon setting 
up the Klein-Gordon lattices, their breather solutions and associated lin¬ 
earization (consisting of both internal modes and wave continuum), we 
explore the notion of their Krein signature. Subsequently, we use asymp¬ 
totic expansions to reveal the nonlinear resonant mechanism leading to 
the main instability result of the present work. In section 3, we numer¬ 
ically corroborate the analytical results. We report computations of the 
Floquet spectrum associated with the breathers in two examples of the 
Klein-Gordon lattice associated with the Morse and hard 0 4 potentials. 
We also present direct numerical simulations highlighting the outcomes of 
the unstable or stable dynamics of two-site breathers. Finally, in section 
4 we summarize our results and present some directions for future work. 


2. Mathematical formalism 

We start with the one-dimensional Klein-Gordon (KG) lattice equation: 

fin + U ( Up ,) — 1 T Un— l); ^ (1) 

where V : M —>• R is an on-site (substrate) potential and e > 0 is the 
coupling constant. The amplitudes of coupled nonlinear oscillators on 
lattice sites form a sequence {u n } n& z € M z denoted by the vector u, 
which is typically defined in the sequence space £ 2 (Z). The solution u 
of the discrete KG equation (jT]) is a function of time t. A local solution 
u exists in C 1 ((—T,T);f 2 (Z)) for some T > 0 if V' is Lipschitz, thanks 
to the Picard’s contraction method and the boundedness of the discrete 
Laplace operator in the lattice equation (fTj). 

For numerical experiments, we consider two prototypical examples of 
smooth on-site potentials, namely the Morse and <f> 4 potentials, which are 
given, respectively, by 

(i) V(u) = ^(e~ u - l) 2 and (ii) V(u) = ^u 2 + ^u 4 . (2) 

We note that K(0) = U'(0) = 0 and U"(0) = 1 due to our normalization. 
The Morse potential (i) in ([2]) is classified as a soft potential because the 
period T of oscillations in the nonlinear oscillator equation 

<P + V'(<P) = 0 (3) 

increases with the oscillator energy E = ^ 2 +U(y?) [I2j. The 0 4 potential 
(ii) in (J2I) is a hard potential because the period T decreases with E. Fig. 


4 


J. Cuevas-Maraver, P.G. Kevrekidis and D.E. Pelinovsky 




Figure 1. Dependence of the period T with respect to the energy E for 
the nonlinear oscillator ([3|) in the Morse (a) and cj) 4 (b) potentials. 


[T| shows the dependence of the period T with respect to the energy E for 
the Morse (a) and <f> 4 (b) potentials. 

The discrete KG equation Jl]) is associated with the Hamiltonian 

H = + ~ u n) 2 + ^ V (u n ). (4) 

nEZ nEZ nEZ 


The energy conservation (j3J) can be used to extend the local solution u € 
C 1 ((—T, T);£ 2 (Z)) to the global solution u € C' 1 (M; £ 2 (Z)) if the potential 
V possesses some coercivity. This is definitely the case for the hard 4> 4 
potential (ii) in (f2j), which supports global solutions of the KG lattice 
equation (JT|) with time-independent bounds on 11u11 ^2 and ||u||f 2 . Due to 
positivity of V for the soft Morse potential (i), the local solution u € 
C 4 ((—T, T);£ 2 (Z)) is extended for e > 0 to the global solution with time- 
uniform bounds on ||u||p and ||<fu||^ 2 , where 5 is the difference operator 
defined by (5u) n = u n+ 1 — u n , n € Z. By using ^||u ||^2 < ||u||^ 2 , we 
obtain a global but linearly growing bound on ||u||^ 2 . 

Let us mention in passing that the soft 1 j) A potential V(u) = \u 2 — \u 4 
supports the blow-up in a finite time for sufficiently large initial data [lj. 
To avoid the finite-time blow-up, we take the Morse potential V in (i) as 
a representative of the class of soft potentials. The Morse potential has 
also applications in physics of DNA eh. 

We consider the breather solutions of the discrete KG equation dTJ, 
which are T-periodic in time and exponentially localized in space. These 
solutions are given by the vector u (t) = u (t + T ) such that u € 
Cper((0> T); £ 2 (%j)) if V is smooth. Breather solutions are constructed 
by the arguments based on the implicit function theorem for sufficiently 
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small e > 0 m (see also review in [29]). In what follows, we always 
assume that the potential V is smooth. 

Adding a perturbation w € C°°(R; £ 2 (Z)) to the breather solution 
u(t) = u (t + T), where u € C^ r ((0,T);£ 2 (Z)), we obtain the linearized 
discrete KG equation in the form 

w n + V"(u n )w n = e(w n+ i - 2w n + w n -i), n € Z. (5) 

According to the Floquet theory, we are looking for solutions of the lin¬ 
earized equation <[5|) in the form w (t) = e A *W(f), where A € C is a spectral 
parameter and W £ C^ r ((0, T)\ £ 2 {Z)) is the eigenvector of the spectral 
problem 

W n + 2XW n + X 2 W n + V"{u n )W n = e(W n+1 -2W n + W n - 1 ), n £ Z. (6) 

Given an admissible value of the spectral parameter A, we can compute 
the Floquet multiplier q by q = e ^ T . The breather u is spectrally stable if 
all admissible values of the Floquet multipliers /i belong to the unit circle. 
Due to the reversibility and reality of the linearized discrete KG equation 
m, the Floquet multipliers /i occur either in pairs on the unit circle {/q q} 
or real pairs {/q/i -1 } or quartets off the unit circle {/q /q /x 1 , q^ 1 }. The 
breathers are spectrally unstable if real pairs or complex quartets exist 
outside the unit circle. 


2.1. Wave continuum 

The wave continuum is defined by solutions of the limiting spectral prob¬ 
lem 

W n + 2XW n + X 2 W n + W n = e{W n+ \ — 2W n + W n _i), n £ Z, (7) 

which corresponds to the zero solution u = 0. Performing the discrete 
Fourier transform, 

/ 7T 

W(t, 6)e ine d0, n £ Z, (8) 

-7r 

and using the Fourier series for the T-periodic functions 

27T 

W(t,0) = Wm(9)e im “°\ u 0 = —, (9) 

mGZ 

we obtain the dispersion equation for the spectral bands 


A = i(±oj(Q) — muo), m £ Z, 


(10) 
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where the function u : [—7r,7r] —» [1, \/l + 4e] represents the fundamental 
band of the wave continuum given by 

w(0) : = ]j l + 4esin 2 6 >g[-7t,7t]. (11) 

In terms of the Floquet multiplier p = e^ T , only two spectral bands show 
up symmetrically on the unit circle. The wave continuum corresponds to 
the Floquet multipliers at p = e ±ia, ( 0 ) T , Q g [—7r, 7 t]. The two spectral 
bands shrink to the two points of infinite multiplicities at /r = e± lT in the 
anti-continuum limit e —» 0. 


2.2. Krein signature 


The Krein quantity represents the symplectic structure of the linearized 
discrete KG equation 0. It is used to characterize Floquet multipliers 
on the unit circle and stability transitions of discrete breathers mm 
Eg. Recently, the same quantity is used in the context of stability of 
the periodic travelling waves in the FPU lattices [7]. Using variables 
{w n ,p n } n £ z, we can write 0 in the form, 


dw n 

dt 


dU dp n 
dp n 1 dt 


on 

dwJ 


n G Z, 


where T~L is the second variation of the Hamiltonian H in 


( 12 ) 

given by 


ft = \^P 2 n+ ^^{ w n+l ~ Wnf + \^V"{U n )w 2 n . (13) 


r) / j Jr TL ' ^ Z. _ J V io~r± uj 1 

nEZ nEZ 


Then, the symplectic (Krein) quantity is given for any w,p£ £ 2 (Z) by 

K = i^2(p n w n -p n w n ). (14) 


nE Z 


We note that K is real and constant in time t. It is identically zero at the 
real-valued solutions of the linearized discrete KG equation (f5|) . However, 
it is nonzero for complex-valued solutions, e.g. for eigenfunctions w (t) = 
e* f2< W(f), where A = iQ and W G C^. r ((0,T); l 2 (Z)) are solutions of the 
spectral problem (|6|). In this case, the Krein quantity K can be written 
in the equivalent form 

K = 2nJ2\W n \ 2 + iJ2(W n W n -W n W n ). (15) 

nE Z nE Z 


For the wave continuum, given by the Fourier transform (jSJ) and the 
Fourier series 0 for H = ±w((9) — muo with 6 G [—7r,7r] and m G Z, we 
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obtain 


/ 7T 

oo{d)\W m (6)\ 2 de. 

-7T 


Since u{9) > 0 for every 9 6 [— 7r, 7r] , according to the dispersion relation 
CD- we have K > 0 for each m-th spectral band with ft = u(6) — mwo 
and K < 0 for each m-th spectral band with 9 = —u(9) — mu o, for every 
m E Z. 


2.3. Internal modes of multi-site breathers 

The internal modes correspond to isolated Floquet multipliers on the unit 
circle that occur commonly in the linearization of multi-site breathers m 
ED. The internal modes and their Krein signature are well approximated 
near the anti-continuum limit of weak coupling between the nonlinear 
oscillators. Asymptotic expansions for eigenvalues and eigenfunctions of 
the spectral problem ([6]) for small e are reported in the most general case 
hr [ 22 ] - Here we review these expansions in the particular case, when the 
nonlinear oscillators are excited on adjacent lattice sites as e —» 0 (see 
also Ham). 

Let us represent the multi-site breather by its limiting configuration 

N 

uiim (t) ■= limii(t) = ^2<p(t + (<Jj - l)T/4)ej, (16) 

3 =1 

where N is the number of excited oscillators, e ? is the unit vector sup¬ 
ported at the j-th lattice node, aj € {+1,-1} for j = 1,2, ...,1V, and 
ip € Cp£ r (0, T) is a unique, even T-periodic solution of the nonlinear oscil¬ 
lator equation ([3]) such that y?(0) > 0. The two nonlinear oscillators at the 
j-th and (j+l)-th lattice nodes are in phase if &jUj +i = 1 and anti-phase if 
(J 3 aj+i = —1. Under this choice of in-phase and anti-phase configurations, 
it is proved that the multi-site breather solution u € C“ r ((0,T);£ 2 (Z)) is 
uniquely continued from its limit (1161) for small e > 0 and the correction 
term u — ui im is Op (e) [20,, 29]. 

The spectral stability problem ([6]) with e = 0 associated with the 
limiting breather (1161) has the Floquet multiplier p = 1 of algebraic mul¬ 
tiplicity 2N and the pair of Floquet multipliers p = e± lT of infinite al¬ 
gebraic multiplicities. The two spectral bands of the wave continuum on 
the unit circle bifurcate from the pair of Floquet multipliers p = e ±iT . 
The Krein signature of the two spectral bands is sign(iF) = ±1. Isolated 
Floquet multipliers may also bifurcate from the pair of Floquet multipli¬ 
ers p = to share the same Krein signature as the one for the spectral 
band they bifurcate from. In the context of the discrete NLS equation, 
the absence of isolated Floquet multipliers bifurcating from the spectral 
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bands can be proved under some restrictive conditions [28]. In either case, 
our main interest is to study internal modes arising as a result of splitting 
of the Floquet multiplier fi = 1 of algebraic multiplicity 21V under the 
perturbation terms with small e. 

Representing the small spectral parameter A = e = A in the spectral 
problem Q, the authors of |29] justified the following asymptotic expan¬ 
sion for the eigenvector W of the spectral problem ([6]): 

N N 

W(t) = £ Cjtp(t+(aj-l)T/4:)ej-£2 Cjv(t+(aj-1)T/4)ej+eW rem (t), 

3 = 1 3 =1 

(17) 

where the even T-periodic function v € C^ r (0, T) is uniquely determined 
by the solution of the inhomogeneous equation 

v + V"(<p)v = 2 ip, (18) 


c = (ci,C2, ...,Cjv) is the vector of projections to be determined, A is the 
new spectral parameter, and W rem is the remainder term of the asymp¬ 
totic expansion. The values of Oj € {+1,-1} for j = 1,2,...,1V are the 
same as in d). 

In the case of symmetric potentials V satisfying V(—u) = V(u), it is 
shown in [29] that W rem £ C' ] ^ r ((0,T);£ 2 (Z)) exists only if the spectral 
parameter A is defined by the matrix eigenvalue problem for the eigen¬ 
vector c: 


T(E) 2 
ST'(E) 


A 2 c = Me, 


(19) 


where S = ip 2 (t)dt is a positive constant coefficient, T'{E) is the 
derivative of the period of the nonlinear oscillator equation @ with re¬ 
spect to its energy E = \dp 2 + V{ip), see Fig. [T] and M is the tri-diagonal 
matrix given by the elements 


Mij — 


a j( a j+i "h &j— l), * Jj 

1, i = j± 1, 


subject to the Dirichlet boundary conditions. In the case of hard poten¬ 
tial with T'(E ) < 0 and for the two-site in-phase breather with N = 2 
and (7 1 = <72, the matrix eigenvalue problem ( 1191 ) has the double zero 
eigenvalue A = 0 and a pair of purely imaginary eigenvalues 

A — ±<Slo. (20) 

The same conclusion holds for the soft potential with T'(E) > 0 but 
applies to the two-site anti-phase breather with N = 2 and oq = —< 72 . 
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The asymptotic approximation (1201) corresponds to the pair of Floquet 
multipliers on the unit circle, which are given at the leading order by 

H = e ±ie ^ n °. 

We shall now substitute the asymptotic expansion (1171) with A = i£lo 
into the Krein quantity m and compute the Krein signature of the 
internal modes (1201) to the leading order in e. After straightforward com¬ 
putations, we obtain K = e* K\ + O(e), where 

N 

K\ = 2flo A\cj\ 2 , A := ip 2 — ipv + v(p. (21) 

i = 1 

By the time-conservation of K, it is clear that A is independent of t. Let 
us compute this quantity explicitly by using the exact expression for v 
satisfying the linear inhomogeneous equation (1181) obtained in [29] : 

v(t) = tip(t) + ^j^d E (p(t] E), 

where <p(t;E) denotes the family of even T(K)-periodic solutions of 
the nonlinear oscillator equation Q satisfying the constraint 0; E) = 
a(E) > 0, where a(E) is the first positive root of V(a(E)) = E. Substi¬ 
tuting v into the formula for A, we obtain 

A = ma E¥ ,(t;E) - ^g E ,p(t E)} = ( 22 ) 

where the last identity holds because the time-independent Wronskian 
of two linearly independent solutions of the homogeneous equation u + 
V"(<p)u = 0 is given explicitly by 

W := <p(t)d E <p(t-, E) - tpd E ip(t; E) = v3(0)«9 E 99(0; E) - ipd E ^(0] E) 

= <f(0)a'(E) = -V'(a(E))a'(E) = -1. (23) 

Substituting (12^1) into (l^ID . we obtain 


K = — 


2n 0 T(E) 
T'{E ) 


£2 


N 

E 

3= 1 


c j | 2 + C , ( e )- 


(24) 


For the eigenvalue A = is^Qo + 0(e) with flo given by (1201) associated 
with the the internal mode of the two-site breather, we have K > 0 for 
the hard potential with T'(E ) < 0 and K < 0 for the soft potential with 
T'(E) < 0. 

Recall that T(E) —>• 27r as E —>• 0, where E —>• 0 represents the small- 
amplitude limit in the nonlinear oscillator equation ([3]). We assume that 
T\E ) remains sign-definite for all admissible values of the energy pa- 
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rameter E. Let us now summarize on the Krein signatures for the wave 
continuum and the internal mode of the two-site breathers: 

• For the hard potential with T'{E) < 0, the internal mode with the 

Floquet multiplier p = e lT ( £2 ^o+Op)) ]~ ias |} ie positive Krein signa¬ 
ture and so is the Krein signature for the wave spectrum associated 
with the spectral band of Floquet multipliers p = e lT ( 1+CJ ( e )). Since 
T(E) < 27r for the hard potentials, we have two situations: 

— 0 < T < 7r: the Krein signatures of the internal mode and the 

wave spectrum in the upper semi-circle coincide; 

— tt < T < 2n: the Krein signatures of the internal mode and 

the wave spectrum in the upper semi-circle are opposite to each 
other. 

• For the soft potential with T'{E ) > 0, the internal mode with 

the Floquet multiplier p = e * T ( £2 ^°+ c, L)) j ias ^he ne g a tive Krein 
signature, which is opposite to the Krein signature for the wave 
spectrum corresponding to the spectral band of Floquet multipliers 
p = e ? ' T ( 1+0 ( e )). Since T(E) > 27r for the soft potentials, we shall 
restrict out attention to the interval of T between 2ir and 47 t (which 
is the 1 : 2 resonance). If 27 r < T < 47 t, we have two situations: 

— 2ir < T <37r: the Krein signatures of the internal mode and 
the wave spectrum in the upper semi-circle are opposite to each 
other; 

— 37 t < T < 47 t: the Krein signatures of the internal mode and 
the wave spectrum in the upper semi-circle coincide. 

We will show that the nonlinear instability occurs when the Krein signa¬ 
tures of the internal mode and the wave spectrum are opposite to each 
other. At the same time, the two-site breathers are nonlinearly stable if 
the two signatures coincide. 

2-4- Asymptotic expansions 

We describe now the main result of this work. We show with the use of 
formal asymptotic expansions that the multi-site breathers are unstable 
in the discrete KG equation ([Tj) if the Krein signatures of the internal 
mode and the wave spectrum are opposite to each other. To study the 
nonlinear dynamics of the internal mode, we adopt the asymptotic expan¬ 
sions obtained earlier for the continuous and discrete NLS models pa eh. 
We assume the following spectral properties: 

P1 There exists a unique internal mode with eigenfrequency such that 
the Floquet multiplier p = e im is isolated from the spectral bands 


Nonlinear instabilities of multi-site breathers in Klein-Gordon lattices 


11 


located at n = e ±luJ ( 0)T , 6 G [—7r, 7r], where u>(6) is given by (fill) . To 
be mathematically precise, we assume that 

0 G (0, min{l — ko^o, mocoo — Vl + 4e}) , (25) 

where uq = ^ is the breather frequency, ko € N is the maximal 
integer such that 1 — kouo > 0, whereas m o G N is the minimal 
integer such that tuquiq — \/l + 4e > 0. 

P2 The double frequency of the internal mode eigenfrequency belongs 
to the spectral band of the wave spectrum, that is, 

either 2f l G (1—koLOo, y/1 + As—kouo) or 2Q G (mocuo—VT+~4e, mQUjQ— 1), 

(26) 

where wo, ko, and mo are the same as for (|25|) . 

P3 There exists a unique eigenvector u G C“ r ((0, T); £ 2 (Z)) of the spec¬ 
tral problem © with A = 0, where u G C^ r ((0, T); £ 2 (Z)) is the 
breather of the discrete KG equation (P). 

P4 The T-periodic functions u n and V'"(u n )\W n \ 2 are respectively odd 
and even in the time variable t for every n G Z, where W G 
Cper((0> T); ^ 2 (Z)) is the internal mode for the eigenvalue A = iQ 
found from the spectral problem ([6]). 

We introduce a small parameter 5 that stands for the amplitude of 
the internal mode. The small parameter 5 is unrelated with the possibly 
small parameter e. Then, we look for an asymptotic expansion 

U(t) = u(i) + W 1 ) (i) + (i) + <5V 3 ) (t) H-, (27) 

where u G C'-((0,T);^ 2 (Z)) is the underlying breather. We choose the 
first-order correction term in the form, 

u W(t) = c(r)W(t)e int + c(T)W(t)e~ int , r = d 2 t, (28) 

where c is the slowly varying amplitude, W G C-((0,T);£ 2 (Z)) is the 
eigenvector of the spectral problem ([6|) for the eigenvalue A = iQ isolated 
from the wave spectrum, see property (PI). 

Separating the variables for the second-order correction term, we rep¬ 
resent 

u (2 >(t) = c(r) 2 P(f)e 2iat + |c(r)| 2 Q(t) + c{T) 2 P(t)e~ 2int , (29) 

where P,Q G C“ r ((0,T);£°°(Z)) are solutions of the linear inhomoge¬ 
neous equations 

P n +4iQP n —4£l 2 P n +V"(u n )P n = s(P n+1 -2P n +P n . 1 )~V'"(u n )W 2 , neZ 

(30) 
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and 

Qn + V"(u n )Q n = E(Q n+1 -2Q n + Q n _ 1 )-V"'(u n )\W n \ 2 , n G Z. (31) 

Although the linear operator in the linear equation m is not in¬ 
vertible because u G C“ r ((0,T);^ 2 (Z)) is a solution of the homogeneous 
equation, see property (P3), the Fredholm solvability condition is satisfied 

V / V"\u n )\W n \ 2 Undt = 0, 

because u n is odd T-periodic and V'"{u n )\W n \ 2 is even T-periodic func¬ 
tions with respect to the time variable t for every n G Z, see prop¬ 
erty (P4). Consequently, there is a unique even T-periodic solution 
Q G C^ r ((0, T); £ 2 (Z)) of the linear inhomogeneous equation (13T1) . 

The linear operator in the linear equation (13011 is not invertible be¬ 
cause 2O belongs to the spectral band of the wave spectrum, see prop¬ 
erty (P2). Consequently, the bounded even T-periodic solution P G 
C , ~((0,T);*°°(Z)) is not uniquely defined unless the boundary condi¬ 
tions are added as n —> ±oo. To obtain a unique bounded solution for P, 
we specify the Sommerfeld radiation boundary conditions as follows. 

Let (7 = 1 if T G (0,7r)mod(27r) and a = — 1 if T G 

(-7T, 27r)mod(27r). Then, we require that the even T-periodic solution 
P G C'-((0,T);£ oo (Z)) of the linear inhomogeneous equation (f30j) satis¬ 
fies the boundary conditions 

R±(t)= Ihn P n (t)e ±M , (32) 

n —>-±oo 

where R± G C™ r (0,T) are uniquely defined and 9q G (0,7t) is uniquely 
found from the solution of the transcendental equations: 

£7 = 1: u](0q) — kouo = 2Q or cr = — 1 : mo^o — w(@o) = 2S7, (33) 

where co(0) is given by the dispersion relation (flTT) . uo = ^ is the breather 
frequency, ko G N is the maximal integer such that 1 — koUJo > 0, whereas 
mo G N is the minimal integer such that mocoo — V1 + 4e > 0. The 
existence of a unique 8q G (0,7t) follows from the assumption (P2). We 
assume (and this will be a subject of the forthcoming work) that there 
exists a unique even P G C“ r ((0,T);£°°(Z)) satisfying (1301) and (|32l) . 

The Sommerfeld radiation boundary conditions (|32j) have the following 
physical meaning. If cr = 1 and e is sufficiently small, the spectral band for 
the wave spectrum corresponding to the Floquet multipliers [i = , 

6 G [—7r, 7r] is located in the upper semi-circle. Then, as n G ±oo, 
the asymptotic solution (|27j) with the radiation boundary conditions ()32l) 
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becomes 


U n (t) ~ 5 2 ^ ± (t) c 2 ( r ) e ^o+i(a;(0 o )-fco^)t + j R ± ( i ^2( r ) e ±m0o-i(^(0o)-fco^o)A +O ( (5 3 


Recall that do £ (0,7r) and u/(#o) > 0, where the derivative of oj in 6 
determines the group velocity of the linear wave packets. Therefore, the 
linear wave packets far from the breather held propagate outwards from 
n = 0, where the breather is localized. On the other hand, if a = — 1, 
the spectral band for the wave spectrum corresponding to the Floquet 
multipliers /i = e iw ^ T , 6 € [ —7r, 7t] is now located in the lower semi¬ 
circle, so that the Floquet multiplier [i = e 2 *^ for the double frequency 
212 is in resonance with the spectral band corresponding to the Floquet 
multipliers /r = e - luJ ( s ) T : Q g [—7r, 7 t]. Then, as n —» ±oo, the asymptotic 
solution (1271) with radiation boundary conditions (|32D becomes 


U n [t) ~ 5 2 ^R ± it)c 2 iT)e ±ineo+ ^ moU> °-^ Go))t + R ± it)c 2 iT)e^ ineo -^ moUo -^ eo))t )+OiS 3 ). 


After the change in the sign of the Sommerfeld boundary conditions (1 3 21) , 
the linear wave packets far from the breather field still propagate outwards 
from n = 0, where the breather is localized. 

Finally, separating the variables for the third-order correction term, 
we represent 

u ®(f) = c(r) 3 G(t)e 3 *™ + F(t,r)e int + F (t,r)e~ int + c(r) 3 G(f)e _MM , 

(34) 

where F(-,r), G £ G“ r ((0, T); £°°(Z)) are solutions of the linear inhomo¬ 
geneous equations for every r. We only write the problem for F, because 
it yields a solvability condition on the slowly varying amplitude c of the 
internal mode. The function F is found from the linear inhomogeneous 
equation: 


F n +2iQF n —F n -\-V"iu n )F n — eiF n j r i—2F n +F n -i)+F[ n , n £ Z (35) 


where the dots denote derivatives with respect to t and the source term 
H is given by 

H n = -2c (w n + inW n ^-\c\ 2 cV'"iu n )iW n Q n +W n P n )-^\c\ 2 cV m iu n )\W n \ 2 W n , n € Z. 

(36) 

The linear operator in the linear equation (|35D is not invertible because 
W £ C^ r ((0,T);£ 2 (Z)) is a homogeneous solution, see property (PI). 

Therefore, a solution exists for F(-,r) £ Cp§ r ((0, T); £ 2 (Z)) if and only if 
the source term H satisfies the Fredholm solvability condition 


^ / W n (t)H n (t, t) dt = 0, 
1 nez-'o 
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which yields the amplitude equation 

iK-j- + /3|c| 2 c = 0, 

(IT 

where K is the same as the time-independent Krein quantity (1151) . 


(37) 


2 i 


r T 


k = ~yJ 1 W n (W n + iVLW n J dt = 2Qj^ \W n \ 2 +i J2(W»W n -W n W n ), 

nGZ ^ nGZ n£Z 


whereas /3 is the coefficient of the cubic term given by 


n€Z U 


V'"(u n )(\W n \ 2 Q n + W 2 P n ) + \v""{u n )\wX ) dt. 


The coefficient /3 is complex because the bounded but non-decaying vector 
P is complex-valued, according to the radiation boundary conditions (l32l) . 
We are only interested in the imaginary part of f3, which can be computed 
by using the linear inhomogeneous equation (1301) and integration by parts: 


2zIm(/3) 


T E V ' "'K) ( W n P n - W 2 P n ) dt 

1 nez Jo 



Pn (Pn + 4iflP n 



^ E / T [ P n( AP )n - P n(AP) n \ dt, 



4 iQP n 


e(AP) n ) 


dt 


where we have used notation (A P) n = P n+ \ — 2 P n + P n -\ for the dis¬ 
crete Laplacian operator. Note that the integration by parts in t yields 
a vanishing result because P € C“ r ((0, T); £°°(Z)). However, because P n 
does not decay to zero as \n\ —>• oo, we have a nonzero result for Im(/3). 
Indeed, we can write 


P n (AP) n -P n (AP) n 


S n -S n - 1 , S n := P n (P n+1 -P n )-P n (P n+1 -P n ), n <E Z, 


interchange the integration and summation, and use the telescopic sum¬ 
mations to obtain 


21Im(/3) 


2e 

Y 


lim S n 

n—>-+00 


lim S n 

n — y —oo 


dt. 


Substituting the Sommerfeld boundary conditions (132li . we arrive to the 
sign-definite expression, 


2zIm(/3) 


4ea sin(6» 0 ) — 



(\R + (t)\ 2 + \R-(t)\ 2 )dt, 
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where sin(0o) > 0 for 6 q € (0,7r). It follows that Im(/3) / 0, if the 
radiation amplitudes R± in the boundary conditions (I ■'121) are nonzero. 
Multiplying (j37j) by c and subtracting the complex conjugate equation, 
we obtain the rate of change for the squared amplitude: 

= — 4 £crsin(0 o )|c| 4 ^ J (\R+(t)\ 2 + \R-(t)\ 2 ) dt. (38) 

It follows from this equation that the squared amplitude |c| 2 decays to 
zero if sign(K) = sign(cr) and grows if sign(J\) = —sign(cr). This is the 
main result of the asymptotic theory. In the former case, we can anticipate 
that the multi-site breathers are stable in the nonlinear dynamics of the 
discrete KG equation. In the latter case, we can predict that the multi-site 
breathers are nonlinearly unstable, in spite of their linearized stability. 

In the approximation of the small-amplitude breathers with the dis¬ 
crete NLS equations [25], T is close to 2ir, so both hard and soft potentials 
feature nonlinear instability of the two-site breathers. This phenomenon 
was discovered recently for the discrete NLS equation in na, following 
upon the earlier abstract analysis of the continuous NLS equation in mi- 
Note that the conclusion changes drastically when the breather period is 
either smaller than n or bigger than 37r. The former case is observed for 
small-period breathers in the hard potentials. The latter case is observed 
for breathers in the soft potentials near the 1 : 2 resonance. In either 
case, the multi-site breathers are expected to be nonlinearly stable. 


3. Numerical results 

Having explored the fundamentals of the nonlinear instability of multi¬ 
site breathers, we now turn to a numerical examination of the discrete 
KG equation (P) for two potentials in (|2|) , namely the soft (Morse) and 
the hard (^> 4 ) potentials. We choose two-site breathers with the excited 
sites at n = 0 and n = 1 oscillating in anti-phase for soft potentials, 
so that uo(0) = —rti(O), and oscillating in phase for hard potentials, so 
that «o(0) = ui(0). According to section 2.3, these two-site breathers 
possess an internal mode with the small eigenfrequency id = + 0(e), 

where ido is given by (1201) . The corresponding pair of Floquet multipliers 
/i = is located near /r = 1 for small e. The Krein signature, 

sign(AT), of the internal mode for the Floquet multiplier /i = in the 
upper half-circle is +1 if the potential V is soft and — 1 if V is hard. 

In general, independently of the expression for the potential V, the pair 
of phonon arcs on the unit circle (which represent Floquet multipliers for 
the wave continuum bands) does not overlap for small e; one arc is located 



16 


J. Cuevas-Maraver, P.G. Kevrekidis and D.E. Pelinovsky 


Morse potential: loq = 0.85 


Morse potential: wq = 0.65 




Figure 2. Floquet argument 9 restricted on [0,7r] versus e for the two- 
site breathers oscillating in anti-phase for the Morse potential (a,b) and 
in phase for the 4> 4 potential (c,d). Red (blue) lines represent spectral 
modes with positive (negative) Krein signature. 


in the upper half-circle and the other arc is located in the lower half-circle. 
The Krein signatures of the two distinct phonon arcs are opposite to each 
other. The arc in the upper half-circle has sign(K) = — 1 if ir < T < 2tt or 
3 vr < T < 47 T and sign(i^) = 1 if 0 < T < tt or 2tt < T < 3tt, where T is 
the breather period. Recall that T < 2tt for hard potentials and T > 2ir 
for soft potentials. For an earlier discussion of these spectral features, see 
e.g. [ZD[22]. 

If 2fi lies inside the wave continuum, the breathers are subject to the 
nonlinear instability mechanism presented herein, as long as the Krein 
signature of the phonon arc at the upper half-circle is the opposite of the 
Krein signature of the internal mode residing in this half-circle. This im¬ 
plies that for the configurations above, nonlinear instability can arise for 
the intervals tt < T <2t r (hard potentials) and 2tt < T < 3n (soft poten¬ 
tials). However, nonlinear instability cannot be observed in the intervals 
0 < T < tt (hard potentials) and 3tt < T < Att (soft potentials). All 
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four cases are represented in Fig. [2], where the dependence of the relevant 
arguments 8 of the Floquet multipliers /r = e %e is displayed versus the 
coupling constant e. For further experiments, we represent the unstable 
and stable cases as follows: 

• Morse potential: unstable case - e = 0.03 for ujq = 0.85 (T ~ 2.357r); 
stable case - e = 0.04 for ujq = 0.65 (T ~ 3 .087r). 

• (j) 4 potential: unstable case - e = 0.3 for luq = 2 (T = vr); stable case 
- e = 1 for wo = 5 (T = 0.47r). 

In order to observe the emergence of the nonlinear instability, breathers 
are perturbed by adding the internal mode (W(0), W(0)), multiplied by 
a relatively small factor 5, to the breather solution (u(0),u(0)). The re¬ 
ported numerical integration results have been obtained by means of the 
4th order explicit and symplectic Runge-Kutta-Nystrom method devel¬ 
oped in 0112]. This scheme preserves the energy up to a 0(10 8 ) factor 
even for the long integration times (> 10 4 ) used here. 

Fig. [3] shows the outcome for dynamics of the two-site anti-phase 
breather in the Morse potential for e = 0.03 for uq = 0.85 (T ~ 2.357r). 
In this case, the internal mode eigenfrequency is ST ~ 0.1002, whereas 
the nearest phonon band shifted by ko = 1 [see (1251) and (1261) ] is located 
in the frequency range [0.15,0.2083]. As a result, the double frequency 
2fl is inside the phonon band. The two-site breather is perturbed by the 
internal mode multiplied by 5 = 0.2. 

Panels (a) and (b) of Fig. [3] include the breather profile with the 
Floquet spectrum and the internal mode (W,W) at t = 0. Panel (c) 
shows the energy density at the central sites n = 0, 1. The energy density 
is defined by 

■ 2 

hn = — h F(rtn) + ~ [(^n — Un+1 ) T (Un ^n— l) ] j (39) 

so that the Hamiltonian in Q can be written as H = ^2 n£ %h n . Panel 
(d) compares the profile of the perturbed breather at f = 0 to the profile 
of the breather near the end of the simulation. 

From Fig. [3fc, we observe the nonlinear instability of the two-site 
breather with frequency coq = 0.85, which manifests by a sudden decay of 
the energy density at the n = 1 site and the growth of the energy density 
at the n = 0 site. A quasi-periodic single-site breather is formed as a 
result of this instability of the two-site breather. 

Fig. 0] shows the outcome for dynamics of the two-site in-phase 
breather in the hard <f> 4 potential for e = 0.3, 5 = 0.2, and uiq = 2 
(T = 7r). For this case, the internal mode eigenfrequency is H ~ 0.4011, 
whereas the nearest phonon band shifted by mo = 1 [see (125|) and (f26l) ] is 
located in the frequency range [0.56,2]. Again, the double frequency 2H 
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Figure 3. The two-site anti-phase breather, predicted to be nonlinearly 
unstable, in the Morse potential for e = 0.03 for ujq = 0.85: the stationary 
breather profile and Floquet spectrum (a), the internal mode (b), the 
energy density versus time (c), and the profile of the perturbed breather 
at t = 0 and at the end of the simulation (d). 
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Figure 4. Similar to Fig. [3] but now for the two-site in-phase breather 
in the hard 4> A potential for e = 0.3 and uq = 2. 


is inside the phonon band. As a result, similarly to Fig. [3] we observe 
the nonlinear instability of the two-site breather and its transformation 
to a quasi-periodic single-site breather. 

We point out that the final state in both computations of Figs. [3] 
and [4] consists of quasi-periodic oscillations of a single-site breather. This 
behavior appears to be generic for the range [ 7 r, 27 t] of the periods T 
supported by the hard (/> 4 potential and the range [27r, 37 t] of the periods 
T supported by the soft Morse potential. 

Figs. [5] and [6] show the cases predicted to be nonlinearly stable in the 
dynamics of two-site breathers. For the soft Morse potential, we take 
the anti-phase breather for s = 0.04 and wo = 0.65 (T « 3.087 t), with 
the perturbation parameter 5 = 0.025. The internal mode possesses the 
frequency fi = 0.1319 and the nearest phonon band is shifted by m 0 = 2 
to [0.22,0.3], so that 2is inside the phonon band. For the hard 0 4 
potential, we take the in-phase breather for e = 1, 8 = 0.2, and ojq = 5 
(T = 0.47r). The internal mode frequency is 14 = 0.8013 and the nearest 
phonon band is not shifted (fco = 0) and is located at [1, 2.23], so that 2 
is again inside the phonon band. In both cases, the internal mode and 
the phonon band in the upper half-circle have the same Krein signature. 
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Despite performing long-time dynamical simulations, an instability of 
the two-site breathers does not seem to arise, confirming our theoretical 
predictions. Due to the Hamiltonian nature of the dynamics and the ex¬ 
citation of the initial perturbation, we observe a quasi-periodic dynamics 
of the original two-site breathers, however, there are fundamental differ¬ 
ences observed between the oscillatory dynamics of Figs. [5] and [6] and 
the unstable ones involving growth (and the eventual destruction of the 
two-site breather states) of Figs. [3] and SI 


4. Conclusions and Future Challenges 

In the present work, we explored a broad question regarding the stabil¬ 
ity of multi-site breathers in Klein-Gordon lattices. This question arises 
in a wide array of settings where discrete breathers emerge, namely in 
nonlinear dynamical lattices. Following the earlier work on discrete (and 
continuous) nonlinear Schrodinger models [111 I15j . we systematically de¬ 
veloped the notion of the Krein signature for the wave continuum and the 
internal modes associated with the excited breather states. Subsequently, 
we produced asymptotic expansions that illustrated the occurrence of the 
nonlinear instability of multi-breathers, when (second) harmonics of in¬ 
ternal modes resonate with the wave continuum of the opposite Krein 
signature. 

However, we also revealed that contrary to the case of discrete (and 
continuous) nonlinear Schrodinger models, depending on the period of the 
breather, it is possible to identify parametric regimes for both soft and 
hard nonlinearities, where the multi-site breathers are nonlinearly stable. 
We have confirmed all of the above conclusions through direct numerical 
simulations of two on-site potentials, namely the soft Morse and the hard 
i j) 4 potentials. 

This instability poses a number of interesting questions for future stud¬ 
ies. On the one hand, it is relevant to explore theoretically whether gener¬ 
alizations of this mechanism can be numerically observed for higher-than- 
second harmonics (with a suitably modified growth law). Furthermore, 
the mechanism is expected to persist not only for multi-site breather con¬ 
figurations bearing a larger number of sites, but also in higher dimensional 
setups and for structurally more complex excited state configurations, 
such as vortex breathers. On the other hand, the nonlinear nature of the 
instability and its weaker-than-exponential growth poses a substantial 
challenge towards its experimental realization. Finally, a deeper under¬ 
standing of the evolution stages of the instability beyond its onset would 
be an important (albeit potentially formidable) task for the dynamics of 
the nonlinear models of interest. 
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Figure 5. The two-site anti-phase breather, predicted to be nonlinearly 
stable, in the Morse potential for e = 0.04 for wo = 0.65: the stationary 
breather profile and Floquet spectrum (a), the internal mode (b), and the 
energy density of the central site versus time (c). Panel (d) shows the 
evolution of the central sites n = 0 and n = 1. 
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Figure 6. Similar to Fig. [3] but now for the two-site in-phase breather 
in the hard c/> 4 potential for e = 1 for ujq = 5. 
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